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In this paper the kernel for the spectral correlation functions of the invariant chiral random matrix 
ensembles with real (/3 = 1) and quaternion real (/3 = 4) matrix elements is expressed in terms of the 
kernel of the corresponding complex Hermitean random matrix ensembles ((3 — 2). Such identities 
are exact in case of a Gaussian probability distribution and, under certain smoothness assumptions, 
they are shown to be valid asymptotically for an arbitrary finite polynomial potential. They are 
proved by means of a construction proposed by Brezin and Neuberger. Universal behavior at the 
hard edge of the spectrum for all three chiral ensembles then follows from microscopic universality 
for (3 — 2 as shown by Akemann, Damgaard, Magnea and Nishigaki. 
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Since its introduction in nuclear physics jjj, Random 
Matrix Theory (RMT) has been applied successfully to 
many different branches of physics ranging from atomic 
physics to quantum gravity (for a recent comprehensive 
review we refer to Hj). One important common ingredi- 
ent is that eigenvalue correlations appear to be insensitive 
to the details of the underlying Hamiltonian. The success 
of RMT is based on this type of universality, and it is no 
surprise that it has received a great deal of attention in 
recent literature [jU^H What has been 

shown is that spectral correlators on the scale of the av- 
erage eigenvalue spacing are insensitive to the details of 
the probability distribution of the matrix elements. Be- 
cause of its mathematical simplicity most studies were 
performed for complex (3 = 2) Hermitean RMT's. How- 
ever, in the case of the classical RMT's it was shown that 
universality extends to real (3 = 1) and quaternion real 
(P = 4) matrix ensembles QJ^Q- This suggests that re- 
lations between correlation functions for different values 
of f3 which can be derived for a Gaussian probability dis- 
tribution |^8|,[29| might be valid for a wide class of prob- 
ability distributions. The main goal of this paper is to 
establish such general relations. As a consequence, uni- 
versality for the much simpler complex ensembles implies 
universality for the real and quaternion real ensembles. 

In this letter we address the question of microscopic 
universality for the chiral ensembles. These ensembles 
are relevant for the description of spectral correlations of 
the QCD Dirac operator. They also appear in theory of 
universal conductance fluctuations in mcsoscopic systems 
|50|,[n]]. In particular, they can be applied to spectral 
correlations near A = 0. According to the Banks-Casher 
formula |}2j , this part of the spectrum is directly related 
to the order parameter E of the chiral phase transition 
(E = Iim7rp(0)/V, where V is the volume of space time 
and p(A) = J2 k £(A— Afc)). It is therefore natural to intro- 
duce the microscopic limit where the variable u = AUE 
is kept fixed for V — > oo. For example, the microscopic 
spectral density is defined by |33| 



Ps(u) = lim (p( — —)), 



(1) 



where the average is over the distribution of the matrix 
elements of the Dirac operator. Successful applications 
of the chiral ensembles to lattice QCD spectra can be 
found in 

The chiral random matrix ensembles for Nf massless 
quarks in the sector of topological charge v are defined 
by the partition function ]3^,|37| 
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where W is a 2n x (2n+v) matrix. As is the case in QCD, 
we assume that v does not exceed y/2n. The parameter 
2n is identified as the dimensionless volume of space time. 
The matrix elements of W are either real (8 = 1, chiral 
Orthogonal Ensemble (chOE)), complex ((3 — 2, chiral 
Unitary Ensemble (chUE)), or quaternion real (3 = 4, 
chiral Symplectic Ensemble (chSE)). For technical rea- 
sons we only consider finite polynomials potentials V(x). 
The simplest case is the Gaussian case with V(x) = E 2 a; 
(also known as the Laguerre ensemble). 

It was shown by Akemann et al. that, for 3 = 2, the 
microscopic spectral density and the microscopic spectral 
correlators do not depend on the potential V(x) and is 
given by the result [Bq] for the Laguerre ensemble. For 
3 = 2 all spectral correlators can be obtained from an or- 
thogonal polynomial kernel corresponding to the proba- 
bility distribution. In their proof the Christoffcl-Darboux 
formula is used to express this kernel in terms of large 
order polynomials. Microscopic universality then follows 
from the asymptotics of orthogonal polynomials. As a re- 
markable achievement, they were able to generalize the 
relation for Laguerre polynomials 
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a L a n (-) = x- a ' 2 J a (2^) 



(3) 



to orthogonal polynomials corresponding to an arbitrary 
polynomial potential. However, their work cannot easily 
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be generalized to (3 = 1 and (3 — 4. The main result of 
the present work is a relation between the kernels for the 
correlation functions of the chOE and chSE and the ker- 
nel of the chUE. This relation is exact for the Gaussian 
ensembles and is valid asymptotically for an arbitrary 
polynomial potential. For (3 = 4, this relation shows 
universality of the microscopic spectral density and cor- 
relators (for (3 = 1 only a partial proof was obtained). 

The partition function (Q) is invariant under W — > 
U^WV where the matrices U and V with dimensions 
determined by W are orthogonal for (3 — 1, unitary for 
(3 = 2, and symplectic for (3 — 4. This invariance makes 
it possible to express the partition function (^|) in terms 
of the eigenvalues x k of WW* as 

Z N f ,= !j{dx kX l^{ Xi )fe- n ^ k y^), ( 4 ) 
J k 

where the Vandermonde determinant is defined by 
AO,) = U k <i( x k ~ xi) and a = N f - 1 + (3v/2 + (3/2. 

For (3 = 2, the spectral correlation functions can be 
evaluated fl40| by expressing the Vandermonde determi- 
nant in terms of the orthogonal polynomials defined by 



(5) 



where we have introduced the potential <j) a {x) = nV(x) — 
a log ir. By using orthogonality relations it can be shown 
that all spectral correlation functions can be expressed 
in terms of the kernel 



2n-l 



(6) 
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The spectral density is given by K%®(x, x) ex-p(—2<fi a (x)). 
Microscopic universality then follows from the following 
generalization of (H) p3[ 



Inn Jhl*ql a (^,n) 
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J2a(u(t)x) 

(u(t)x/2) 2a '' 
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in the normalization qf, a (0) \J h 2 "- = 1. The function u(t) 
follows from the asymptotic properties of the leading or- 
der coefficients of the q 2 . a {x) and the normalizations /i| a . 
Its value at t = 1 is given by u(l) = 2np(0). 

In order to perform the integrations by means of or- 
thogonality relations for (3 = 1 and (3 = 4, one has to 
introduce the skew-orthogonal polynomials p9|p0| . Be- 
low, we first discuss the case (3 = 1 and then give general 
outlines for the case (3 = 4. 

For (3=1, the skew orthogonal polynomials of the 
second kind are defined by 



< Ri, Rj >r— Jij- 
with the skew orthogonal scalar product 



(8) 



<f,9> R = dxe- 2 ^f(x)Zg(x), (9) 
Jo 

and nonzero matrix elements of Jy given by J2fc,2fc+i = 
— ^2fc+i.2fc = — 1- The operator Z is defined by 



Zg(x) 



dye^e(x-y)e-^g(y), 



(10) 



Here, e(x) = x/2\x\. It can be shown that all correlation 
functions can be expressed in terms of the kernel [p9]j4C|] 



K 1 (x,y) = / dze-^k 1 (z,y)e~ M y\ 



o 

where we have introduced the pre-kernel 

2n-l 

k i(y,z) = ^2 Ri{y)JijRj{z). 

i,j=0 

In particular, the spectral density is given by 
p(x) = K 1 (x,x) - ii^i(oo,a:). 



(11) 



(12) 
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A general scheme for the construction of skew- 
orthogonal po lynomials was introduced by Brezin and 
Neuberger |41j| . The idea is to express them in terms 
of orthogonal polynomials defined by (^. For technical 
reasons we expand in the polynomials q k a+l (with weight 
function x 2a+1 exp(— 2nV(x))). The skew-orthogonal 
polynomials of degree i can thus be expressed as 



R t {x)=Y J T l] qf + \x), 

3=0 



(14) 



where T is a lower triangular matrix with nonvanishing 
diagonal elements. An essential role is played by the 
inverse, L, of the operator X^Z with Z defined in ( |l0| ) 
and Xg(x) = xg{x). It can be easily verified that 



L = X(d-<f>' a (X)) + i. 



(15) 



The matrix representations of the operators X, Xd, 
X~ X Z and L in the basis q\ a+1 will be denoted by Xm, 
Dku Yu and L k i, respectively. In the remainder of this 
derivation the index 2a + 1 will be suppressed. 
In matrix notation (||) can be rewritten as 

TYT T = - J. (16) 

By using that LY = 1, this relation can be expressed as 

L = T T JT. (17) 

It can be shown that the matrix L k i is a band matrix 
with width determined by the order of the polynomial 
potential V(x). It then follows that T is a band matrix 
as well EJ. For example, for a Gaussian potential we 
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have that T 2m ,k = a S 2m ,k and T 2m +i,k = boS 2 m+hk + 
biS2m,k + &2^2m-i,fc with coefScients derived in 

It turns out that we do not need explicit expressions 



for the Tij . The pre-kernel (112) can be expressed as 



(18) 



i,j—0 k<i l<j 



In this relation the indices i and j run up to 2n — 1 in 
contradistinction to the relations ( |l7| ) where they run up 
to oo. However, it follows from the band structure of 
L that the number of terms outside the range in ( O ) is 
of the same order as the degree of the polynomial po- 
tential which is finite. These terms are negligible in the 
continuum limit of the type (Q) where the q 2a+1 (x) and 
Lki depend smoothly on k and / (notice that L k i is not 
smooth in \k — l\). However, for x around zero and y 
near the largest zero of qi{y) we expect potentially non- 
negligible contributions. We thus have that 



2n-l 

h(x,y)- X! <lk(x)L k iqi(y). 
k,i=a 



(19) 



By means of a partial integration the matrix elements of 
1 — X(j)' a (X) in L can be expressed in terms of the matrix 
elements of Xd. This results in 



Lki — 
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zdze-^Xm (z)zdq k (z) - q k (z)zd qi (z)) 



-(D kl -D lk ) 



(20) 



The matrix elements of D can be re-expressed as xd x or 
yd y . We finally arrive at a remarkably simple expression 
for k x {x,y), 



h(x,y) ~ -{ydy 



xd x )K^{x,y) 



(21) 



With the help the asymptotic properties of the qi a (which 
are the same as for the Laguerre polynomials) this rela- 
tion can be further simplified to (up to an overall factor 
determined by the average spectral density) 



d x )K%(x,y). 



(22) 



Universality of the unitary kernel K%®(x, y) at the hard 
edge has been well established [^3| for the chiral ensem- 
bles, whereas universality at the soft edge was shown 
in |l^,|r^]. We therefore expect universal behavior of 
k\{x,y) in these domains. 

Let us finally focus on the spectral density. Using (|1J 
and (E2), for a Gaussian potential it can be expressed as 



p(x) 



> 1 

dye-^Mx - y)-(d y - d x )K%(x,y). 

(23) 



In the microscopic limit where n — > oo at fixed z = xn 2 
the factor exp(— 2nV(x)) — » 1 and K 2 n approaches its 
universal limit. However, in one of the terms contribut- 
ing to the integral the microscopic limit and the inte- 
gration cannot be interchanged. It can be shown that 
there is an additional contribution with x near zero and 
y near the edge of the spectrum. Naively taking into ac- 
count this contribution for nongaussian potentials leads 
to a microscopic spectral density that differs from the 
universal expression. Alternatively, we have established 
universality of the microscopic spectral density by means 
of Monte-Carlo simulations. Apparantly, the smoothness 
assumptions in the derivation of (21) are violated in this 
case. The edge contribution resides in the term Ki(x, oo) 
in ( |l3| ) . In the first term contributing to spectral density, 
Ki(x,x), the microscopic limit and the integral can be 
interchanged. This establishes universality of K\{x,x). 

In the case of a Gaussian potential the edge contribu- 
tion can be obtained from the asymptotic expansion of 
the Laguerre polynomials in this region (an expression in 
terms of Airy functions). In a future publication, we hope 
to establish a possible relation with universal behavior of 
the ql a near the edge of the spectrum ■ 

The above analysis carries through for the symplectic 
ensemble. In this case there are no contributions from 
the soft edge and universality of the microscopic spectral 
density can be shown rigorously. For 8 — 4 (with an ad- 
ditional factor 1/2 in the exponent of (Q)), the correlation 
functions can be expressed in terms of the kernel 



k 4 (x,y) = 



2n-l 



i{x)JijQj(y), 



(24) 



This is the central result of this paper. It is valid asymp- 
totically both at the hard and the soft edge of the spec- 
trum where a continuum limit of the orthogonal polyno- 
mials q k a+1 exists. However, as will be argued below, 
the result (^2|) is not valid for x near the hard edge 
and y at the soft edge of the spectrum. This result 
relates the orthogonal pre-kernel to the unitary kernel 
K 2 n(x,y) which has been studied elaborately in the lit- 
erature p3|-p5[ . The relation ( |22] ) is exact for a Gaussian 
potential in which case it coincides with the result ob- 
tained in g,||,|6|jT^g|. 



where the Qi(x) are skew orthogonal polynomials of the 
first kind which are defined by the skew-scalar product 



< f,9 >Q- 



dx 



e-**-Mf{x){L-i)g(x), (25) 



with the operator L defined in (|l5|). In this case we 
express the Qi(x) in terms of the polynomials q 2a ~ 1 (x), 



Q k (x) = Y, S ^ 1{ 



(26) 
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The matrix elements of the operators are also in this 
basis. From the orthogonality relation < Qk,Qi >Q= Jki 
it can be shown that SLS T = — J from which we derive 
S T JS = XYX~ X . Again, due to the band structure of 
Lki, the range of the summations in this relation and 
in ( p4| ) differs by a finite number of terms which can be 
neglected in the continuum limit. We thus find 

2n-l 

k 4 (x,y)c ]T ql a -\x)[zX-') ki q^-\y) 



k,l=0 
oo i 

az 



e(y~z)K%-\x,z). (27) 



o 



z 



Universality of ki(x,y) thus follows from universality of 
K2n~ 1 (x, z). This relation is exact for a Gaussian poten- 
tial and reproduces the result found in 29 1. 

In conclusion, we have shown that relations between 
the kernels for the chOE and chSE and the kernel for 
chUE are not accidental but follow from an intriguing un- 
derlying mathematical structure. Under certain smooth- 
ness assumptions these relations are valid asymptotically 
for an arbitrary polynomial potential. Microscopic uni- 
versality for f3 = 4 and in part for (3 = 1 thus follows 
from universality at /? = 2 at hard edge of the spectrum. 
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